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When a classical device suddenly perturbs a degenerate Fermi gas a semiclassical non-equilibrium 
Fermi state arises. Semiclassical Fermi states are characterized by a Fermi energy or Fermi mo- 
mentum that slowly depends on space or/and time. We show that the Fermi distribution of a 
semiclassical Fermi state has a universal nature. It is described by Airy functions regardless of the 
details of the perturbation. In this letter we also give a general discussion of coherent Fermi states. 



1. Introduction Among various excitations of a de- 
generate Fermi gas, coherent Fermi states play a special 
role. A typical (not coherent) excitation of a Fermi gas 
consists of a finite number of holes below the Fermi level 
and a finite number of particles above it. Instead, a co- 
herent Fermi state involves an infinite superposition of 
particles and holes arranged in such a manner that one 
can still think in terms of a Fermi sea (with no holes in 
it), the Fermi level of which depends on time and space. 

Coherent states appear in numerous recent propos- 
als about generating coherent quantum states in nano- 
electronic devices and fermionic cooled atomic systems. 
These states can be used to transmit quantum informa- 
tion, test properties of electronic systems and to generate 
many-particles entangled states. 

Coherent Fermi states can be obtained by different 
means. One is a sudden perturbation of the Fermi gas. 
For example, a smooth potential well, the spatial extent 
of which much larger than the Fermi length, is applied to 
a Fermi gas. Fermions are trapped in the well. Then the 
well is suddenly removed. An excited state of the Fermi 
gas obtained in this way is a coherent Fermi state. This 
kind of perturbation is typical for various manipulations 
with cooled fermionic atomic gases. 

A realistic way to generate coherent Fermi states in 
electronic systems is by applying a time dependent volt- 
age through a point contact, typifying many manipula- 
tions with nanoelectronic devices [TJ [2] . 

Although the realization of coherent states in elec- 
tronic systems experimentally is more challenging than 
in atomic systems, we will routinely talk about electrons. 

From a theoretical standpoint coherent Fermi states 
are an important concept revealing fundamental proper- 
ties of Fermi statistics. Coherent states appeared in other 
disciplines not directly related to electronic physics. Ran- 
dom matrix theory (RMT) [3] , non-linear waves [4] , crys- 
tal growth [5j [6] , various determinental stochastic pro- 
cesses [7 , asymmetric diffusion processes [8], to name 
but a few. 

Unless a special effort is made (see e.g., [2 ) coherent 
Fermi states involve many electrons and are such that 
space-time gradients of the electronic density are much 
smaller than the Fermi scale. These state arise as a result 
of perturbing tye Fermi gas by a classical device. We 



call them semiclassical Fermi states. They are the main 
object of this paper. We will show that semiclassical 
Fermi states show great degree of universality as well as 
their single electron counterpart studied in [TJ |2] . 

A general coherent Fermi state is a unitary trans- 
formation of the ground state |0) of a Fermi gas: 
\U) = U\0), U = e ifz(x) P (x)dx-ifn(x)v(x)dx^ where p 

and v are operators of electronic density and velocity 
[p(x), v(y)] = — — V5(x — y) and and IL(x) are two 
real functions characterizing the state [IJ[9]). 

For the purpose of this paper it is sufficient to assume 
that the motion of electrons is one dimensional and chiral 
(electrons move to the right), although most of the results 
we discuss are not limited to one-dimensional electronic 
gases. To this end, edge states in the Integer Quantum 
Hall effect may serve as a prototype. In a chiral sector the 
operators of density and velocity are identical v\ right) = 
right). The contribution of the right sector is easy to 
take into acoount. In this case the Fermi coherent state 

\U) =e */*(*M*)^| > (l) 

is characterized by a single function <1>. 

The function <1> can be undersstood as the action of 
an instantaneous perturbation by a potential frupV<I>(:r), 
or, if we ignore the electronic dispersion the function 
is the action of a time-dependent gate voltage eV(t) = 
—H^^(xo-VFt) applied through a point contact (located 
at xo) [16] . 

Fermi coherent states feature inhomogeneous electric 
density p(x) and current I{x). Their expectation values 
give a meaning to the functions S and II. In the chiral 
state, where the electric current and the electronic den- 
sity are proportional, their expectation values are gradi- 
ents of <1> 

h(U\p(x)\U) = —(U\I(x)\U) = p F + fiV$. (2) 

Consequently hV$(x) plays the role of the space-time 
modulation of the Fermi point, such that all states with 
momentum (or energy) less than Pf{x) = HV$(x) (or 
Ep(x) — vphVQ) are occupied (no holes in the Fermi 
sea) as shown in Fig. [I] We shall refer to the function 
Pf(x) as the Fermi surface, with some abuse of nomen- 
clature. We shall also term the region in phase space 
around Pf(%) as the 'Fermi surf. 
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FIG. 1: Modulated Fermi Edge. The Wigner function and 
Fermi number are almost 1 below the Edge (shaded area) 
and vanish above the Edge, having universal character in the 
'surf. 




FIG. 3: Universal behavior of the Wigner function vs. mo- 
mentum for a concave Fermi surface obtained from (13]). 




FIG. 2: Universal asymptotes of the Fermi number. The 
graph (blue curve) shows the Fermi number in units of h/£ 
for the example (p(x)) = po + (n / £) cos(x / £) computed from 
Eq.(|17[). Black curves are the asymptotic forms obtained 
from the universal Fermi number formula Q. The univer- 
sal asymptotes are magnified. 



The question we address in this letter is: what is the 
Fermi distribution in the 'Fermi surf of a modulated 
Fermi point, namely around and between extrema of 
P F {x)l 

We will be especially interested in semiclassical Fermi 
states. These states involve many excited electrons. 
Its typical range is larger than the momentum spacing 
W<£ ^> A (but still is smaller than the Fermi momen- 
tum). 

We show that, quite interestingly, the semiclassical 
Fermi surf features a universal Fermi distribution. There 



the Wigner function (11) is described by the function 
Ah(s) = J™ Ai(s')ds': 



uf(x^p + Pf(x)) ~ Aii 



hip 



where the scale k = \H 2 Pp(x*)/2\~ 1 / 3 , and the offset 
Pf{x) are the only information about the state that en- 
ters the formula. This formula holds close to any point 
where the Fermi surf is concave P^(x) < 0. The Fermi 



occupation number (12) also displays universal behavior. 



Near a maxima, but well above the minima of the surf 
the Fermi number reads 

n F (P F (x*) +p) « ftA [[Ai'(/^)] 2 - («p)Ai 2 («p)l . (4) 

If the Fermi surf is convex rather than concave, then par- 
ticle hole symmetry tif(p,x) — >• 1 — 7ip(— p, provides 
the result for the Fermi number and Wigner function. 



Figures (T|2) illustrate the universal regimes. Here A is 
the momentum spacing (27rS/A is the system volume). 

The goal of the letter is twofold: to emphasize these 
simple, albeit universal, distributions, and, also to collect 
a few major facts about Fermi coherent states. 

2. Coherent Fermi states The formal definition of 
a Fermi coherent states starts with the current algebra 
(see e.g., [11 ). To simplify the discussion and formulas 
we consider only one chiral (right) part of the current 
algebra. 

Current modes are Fourier harmonics of the electronic 
density p(x) = ^Zk>o e ^ kx Jk- An electronic current 
mode Jk = J2 p c^Cp+k ( we count electronic momentum 
from the Fermi momentum Pf), creates a superposition 
of particle- hole excitations with momentum k. Positive 
modes annihilate the ground state, |0), a state where 
all momenta below pp are filled: Jfe|0) =0, k > 0. 
Negative modes are Hermitian conjugated to the posi- 
tive modes J-k = (J+fc)^. 

Chiral currents obey a current (or Tomonaga) algebra: 



[Jk, Jl] = -^5k+l,Q- 



(5) 



A Fermi coherent state \U) is denned as an eigenstate of 
positive current modes: 



Jk\U)=Pk\U), k>0. 



(6) 



(3) As follows from (!]), p^ = A J e *> kx 'd^(x) / '(2ir) are pos- 



itive Fourier modes of the function $(x). Assuming that 
the total number of particles in the coherent state is the 
same as in the ground state, or that the dc component 
of current is Iq = ^j^Pf, i.e., J d$> = 0, we obtain 

\U) = Z~ 1/2 e^ k >° ^ J - fc |0), Z = e^ k >° . (7) 

Using normal ordering with respect to the ground state 
(where all positive modes of the current are placed to the 
right of negative modes) the unitary operator reads [T7] : 

e % J <5>(x)p(x)dx _ %-\ :e if $(x)p(x)dx ._ 

= z~ie^ k >° iP kJ - k e ~^ k >o ip~ kJk . (8) 
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A coherent state represents an electronic wave-packet 
which is fully characterized by the electronic density. It is 
a simple exercise involving the algebra of the current op- 
erators to show that the function V3>(x) is a non-uniform 
part of the density as is in Alternatively, one can 
use Z as a generation function (p(x)) = (U\p(x)\U) = 
p + 2Re'£ k>0 ke ikx dp h logZ. 

Coherent states obey the Wick theorem. The Wick 
theorem allows to compute a correlation function of 
any finite number of electronic operators, as a deter- 
minant over the one-fermionic function K(xi,x 2 ) = 
(U\^(x 1 )^(x 2 )\U), where ^{x) = (h/A) 1 / 2 £ p e^%is 
an electronic operator. The one-fermionic function can 
be computed with the help of the formula: 



which leads to the expression: 



-i<&(x) 



ip(x) 



(9) 



e i®(x 1 )-i®(x 2 )-jr PF (x 1 -x 2 ) _ 

K(xi,x 2 ) = 7} v , (10) 

l(X! - X 2 ) 

valid for A\x± — x 2 \ <C H. An equivalent object appears in 
RMT where it is often called - Dyson's kernel. We adopt 
this name. As points merge one recovers the density ([2| 
K(x,x) (p(x)) V$. 

3. Wigner function and Fermi occupation number 
The Wigner function is defined as Wigner transform of 
the Dyson kernel 



n F (x : p) 



x+ y -,x- y -)e-^Uy (11) 



The meaning of the Wigner function is clarified away 
from the surf. There it means an occupation of electrons 
in the phase space (x,p): 1 below a surf, above. On 
the surf Wigner function is not necessarily positive. 
The Fermi number 



A f 

n F (p) = (U\clcp\U) = — / n F (x,p)da 

is the Wigner function averaged over space. 
Combining ( 10 ) and (fTTT) we write 



(12) 



n F (x,p) 



2iri 



y-iO 



(13) 



where we denoted P F {x) = hV<& = (p(x)) as in (|2|. 

Below we evaluate the integral ( fl3| ) semiclassically 
bearing in mind that $ is of a finite order as h —> 0. 

A universal regime arises at the Fermi surf, n\p — 
P F (x)\ c± 1. In this case it is sufficient to expand P F (x) 
in a Taylor series around extrema of P F (x) to second or- 
der P F (x) = P F (x*) + \P£(x*)(x - x*) 2 + . . . . Then the 
integral (13) becomes the Airy integral given in Eq. (|3j. 



Further integration over space yields Q. 

In this regime the Dyson kernel in the momentum 

space K plyP2 = (U\c^ c P2 \U) reads: 



Ai(^pi)Ai / (^p 2 ) - Ai(^p 2 )Ai / (^pi 
(Pi -P2) 



(14) 



This is the celebrated Airy kernel appearing in numerous 
problems as the limiting shape of crystals [5] , asymmetric 
diffusion [8 , edge distribution of eigenvalues of random 
matrices [13], etc. 

The Fermi number Q can be directly obtained from 
the kernel by taking a limit p\ — »• p 2 in (14). At 
large positive momenta {np — >• +oo) the Fermi num- 
ber behaves as n F (p) 



_A_ e -f(*p) J 
8irp 



and 



-pK 



^ cos(|(— npf*! 2 ^ for large negative mo- 



menta within the surf. 

Away from the universal region of the surf the Fermi 
distribution can be computed within a saddle point ap- 



proximation. The saddle point of the integral (13) is: 



P F (x- 



2 } 



Pf(x- 



2p 



(15) 



It has pairs of solutions ±y*(x,p). Let P ma x = 
max(Pp(x)) and P m in = mm(P F (x)) be adjacent ex- 
trema of the surf. Without loss of generality we may 
assume that (x,p) is outside the Fermi sea p > P F (x). 
The particle hole symmetry n F — >• l — n F helps to recover 
the case when the momentum is inside the sea. If p is in 
the surf, p G (P m in, Pmin), then some saddle point pairs of 



( 15 ) may be real. Their contribution produces oscillatory 



features with a suppressed amplitude. If p hovers above 
the surf, p > P max , then the saddle points are imaginary. 
Their contributions are exponentially small. 

Between two adjacent extrema the Wigner function 
reads: 



n F {x,p) 



8n\y* 



2 sin (ft 



pe 



(Pmin i P max) 
P ^ Pmax 5 



where HQ = — Jp F ^ y*{x,p')dp' . In the surf it is half the 
action of a semiclassical periodic orbit - the area of the 
graph y*(x,p) vs p. In the universal regime, when one 
approximates y*(x,p) ~ Hk, 3 / 2 (p — P F (x)) 1 ^ 2 this equa- 
tion reproduces asymptotes of Eq. (J3| : n F (x, p-\-P F (x) ) ~ 

(87r)- 1 /2( Kp )-3/4 e -3(« P ) 3 / 2 _ 

A Fermi coherent state with a periodic current is an 
instructive example. It corresponds to "quantum pump- 
ing" - periodic transfer a charge through the system by 
applying a periodic voltage through a point contact. Set- 
ting (p(x)) = po + (n / £) cos(x / £) the Dyson kernel in 
the momentum representation becomes the integer Bessel 
kernel 

2 mi —m 2 

where p± j2 = tvm\, 2 /£, where mi ^ m 2 are integers. The 
Fermi number is given by (m F = p F £/h): 



n F (p) = - -sign(p) 



J 2 (n) 



(17) 
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This formula allows to compare the asymptotes near 
the edges to the universal expression above. Using 
the homogeneous asymptote of Bessel function J m (m — 
(m/2) 1 / 3 ~ (2/m) 1 / 3 Ai(C) at large m, one recovers 
(14) and Q. Fig. [2] illustrates the universal asymptote. 

4. Holomorphic Fermions as coherent states To con- 
trast semiclassical coherent Fermi states and quantum 
coherent Fermi states, we briefly discuss special coherent 
states known as holomorphic fermions [T4] . 

Holomorphic fermions are defined as a superposition 
of fermionic modes ip(z) = ^Z p e^ pz c p , with a complex 
"coordinate" Imz < 0. 

Holomorphic fermions are coherent states since they 
can be represented as an exponent of a Bose field - dis- 
placement of electrons <p(z) = Xl/c^o \ e ^ kz ^ 01 E] 



(18) 



A function <I> for a string of fermions 

n?=i (^teMCi)) 10} is = nr=i fEf fE§- The 

density (or current) of these states consists of Lorentzian 
peaks, each carrying a unit electronic (positive/negative) 

charge (p(x)) - po = J2i Im _ ^C*)' so tnat tne 

state is a set of single electronic pulses. For a possible 
applications of these states in nano-devices see [TJ [2]. 
As the complex coordinate approaches the real axis 
£ = x — iO a holomorphic fermion operator becomes an 
electronic operators i/)(z) — »• ip(x) as its density becomes 
a delta-function. 

Coherent states formed by a single holomorphic 
fermion carries a unit charge in contrast to semiclassical 
Fermi states. The Wigner function of this state follows 
from (fTol) 



n F (x,p) 



2iri J — z x — I — z 



y-iO' 



where z = X — ia, X is a real coordinate of fermion and 
a is its width (a A <C h) . a is positive for an annihilation 
operator and negative for a creation operator. Evaluating 
this integral we obtain 



/ \ ^/ \ ^/ \ sin2p(x-X) . 
n F (x,p) 6(— p) - Q(-pa)a v - — -e^ pa . 



n F (p) = e(-p)-6(-pa) 



2aA 



(19) 



Noticeable features of this distribution are: the Fermi 
function jumps on the Fermi edge; beyond the Fermi 
edge, the Fermi number and the Wigner functions de- 
cay exponentially, if a > a holomorphic fermion i/)(z) 
acts like an annihilation operator removing particles from 
Fermi edge (vice versa if a < 0); the Wigner function fea- 
tures Friedel's type oscillation with a distance x — X to 
the center of the fermion. 

6. Fermi coherent states and Random Matrix Ensem- 
bles. We complete the letter by a brief discussion of the 
relation to the theory of Random Matrices. 



Consider n electrons from the ground state of the 
Fermi gas in positions xi,...,x n . We write this state 
as (xi,...,x n | = (0|cJ F . . .cl F _ n x/>(xi) . . .^(x n ), where 
we set momentum spacing to 1 for brevity. Let us ask 
for the probability to find this state in the coherent 
state \U). It is |^(xi, . . . , x n )| 2 , where ^(^i, . . . , x n ) = 
(x u ...,x n \U)/(0\U). 

Fermions in coherent states obey the Wick theorem: 
a matrix element of a particle-hole string inserted be- 
tween two (generally different) coherent states is a Slater 
determinant built out of particle- hole matrix element. 
In particular ^(x 1 ,...,x n ) = det (^p F -i(x m )) lm£n is 
a Slater determinant of a single particle matrix element 

®pf( x ) = (°\ c p F ' l P( x )\ u )/( \ u )' We compute it with the 
help of the current algebra and bosonic representation 
dl8|). Up to a normalization 



*p(*) 



(20) 



-ipz - 



This formula features the complex curve Y p (z) 

^2k>o iPk e ^ kz S a useful characteristic of the coherent 
state. The function Y p (z) + ipz is analytic in the upper 
half-plane. Its boundary value on the real axis is 

Y p (x) +ipx=- (V(x) - i$(x)) = i^yZ^e^ kx , (21) 

k>0 

where V(x) and &(x) are real and imaginary part of the 
boundary value of the analytic function. They are con- 
nected by the Hilbert transform. 

For example, a complex curve for a pumping consid- 
ered in Sec. 3 is Y p (z) = —ipz + |fe?. In the case 
of a string of fermions the curve is Y p (z) = — ipz + 



Computing the Slater determinant of ([20]) we obtain 
. . . ,x n ) 



Z-^ 2 e~i S=i( y (^)-^^)A(x), (22) 



where A(x) = det (e^ pxi ) = rii> 7 (e* Xi 

V / p,l<n J V 

is the VanderMonde determinant. 
The normalization factor in ([22 



/n 
JJ \ e i*i _e^pj| e -^) &i (23) 

n>i>j>l i=l 

is the the partition function of eigenvalues of a circular 
unitary n x n-matrix [3 . At the limit of vanishing spac- 
ing one replaces e^ x — » 1 + ihx. In this case coherent 
Fermi state is described by Random Hermitian Matrix 
ensemble. 

If n is large, Eq. (22) can be interpreted as a co- 
ordinate representation of the coherent state. A Fermi 
coherent state may be thought as a Fermi sea filled by 
particles (without holes) with wave functions (20) and 
p = 1, . . . , n. The coordinate representation provides an- 
other avenue to compute matrix elements discussed in 
this paper as a limit n— ^ 00, pp/n ^ 00. Some of them 
have been studied for various reasons in the theory of 
Random Matrix Ensembles (see e.g., [15 for derivation 
of the Dyson kernel). 
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